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An  Analytical  Soil  Thermodynamic  Model  for  the  Diurnal 
Variation  of  Ground  Surface  Temperature 


I.  IMKODI  CHON 

The  determination  of  diurnal  variations  of  ground  surface  temperature  plays  a 
critical  role  in  the  modeling  of  the  atmospheric  planetary  boundary  layer  and,  by 
extension,  in  mesoscale  numerical  weather  prediction  (NWP)  models.  In  theory, 
the  problem  is  a  relatively  simple  one  because,  to  a  very  good  approximation, 
these  variations  may  be  represented  by  a  linear  partial  differential  equation. 
Assuming  homogeneity  and  constancy  in  thermal  characteristics  of  a  soil  slab  that 
is  subject  to  linear  forcing  at  the  bounding  surfaces,  we  can  easily  solve  the  linear 
heat  equation  analytically.  Ground  surface  temperature  may  then  be  obtained  by 
applying  the  solution  at  the  upper  surface  of  the  slab.  Solutions  of  linear  flow  of 
heat  in  a  slab  bounded  by  two  parallel  planes  have  been  given,  for  example,  by 
C  arslaw  and  Jaeger^  for  various  kinds  of  linear  forcing  at  the  bounding  surfaces. 

In  geophysics,  C  ampbelP  has  given  a  time-dependent  solution  of  the  temperature 
distribution  in  an  infinite  soil  layer  that  is  forced  at  the  upper  surface  by  a  peri- 

(Keceived  for  publication  18  December  1986) 

1.  C'arslaw,  H.S.  ,  and  Jaeger,  J.C.  (1959)  Conduction  of  Heat  in  Solids.  Oxford 

University  Press,  London,  510  pp. 

2.  (  ampbell,  Gaylon  S.  (1977)  An  Introduction  to  Environmental  Biophysics. 

Springer-Verlage,  1  59  pp.~~ 
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odic  temperature  wave.  Blackadar  has  provided  a  solution  for  a  soil  model  of  in¬ 
finite  depth,  subject  to  the  boundary  conditions  that  (1)  the  top  of  the  soil  is  in  in¬ 
stantaneous  thermal  equilibrium  and  is  heated  by  fluxes  that  are  represented  by  a 
single-component  cosine  function  in  time;  (2)  the  temperature  at  infinite  depth  ap¬ 
proaches  a  constant  value  T^j. 

At  the  ground  surface,  the  various  flux  terms  of  the  energy  balance  equation 
unfortunately  are  not,  in  general,  linear  functions  of  Tg,  the  temperature  at  the 
ground.  These  terms  include:  short  wave  flux  due  to  solar  radiation,  infrared  flux 
from  the  atmosphere,  thermal  emission  at  the  ground,  sensible  heat  transfer  be¬ 
tween  the  air  and  the  ground,  latent  heat  flux  due  to  evaporation/condensation  and 
rainfall,  and  soil  heat  flux  due  to  the  vertical  thermal  gradient  within  the  soil.  Of 
these,  only  the  first  has  well  defined,  externally  forced  diurnal  and  seasonal  peri¬ 
ods  and  then  only  under  cloudless  skies.  Therefore,  the  difficulty  in  obtaining  an 
anal3^ical  solution  for  the  heat  equation  lies  not  in  the  differential  equation  itself 
but  rather  in  accurately  prescribing  the  external  forcing  and  in  the  nonlinear  feed¬ 
back  nature  of  the  upper  boundary  condition.  For  this  reason,  values  of  T  in 
NWP-type  models  are  usually  determined  by  numerical  methods,  Deardorff  found 

the  so-called  force-restore  method  to  be  the  most  satisfactory.  This  method  is 

5  6 

currently  used  in  many  NWP-type  models.  ' 

In  this  report,  we  present  a  perturbation  method  that  may  be  used  to  obtain 
approximate  Tg  values  analytically  for  a  soil-slab  model.  The  key  is  to  decompose 
Tg  into  a  basic  part  and  a  perturbed  part  and  then  approximate  all  surface  flux 
terms  other  than  those  due  to  solar  heating  as  linear  functions  of  the  perturbed 
temperature.  Once  such  approximations  are  made,  our  problem  then  becomes 
that  of  linear  heat  flow  in  a  finite  soil  slab  subject  to  boundary  conditions  that  are 
linear  functions  of  the  perturbed  temperature.  Such  a  problem  may  be  solved  ana- 
l3^ically.  Tg  is  then  obtained  simply  by  applying  the  solution  at  the  upper  surface 
of  the  slab.  This  model,  which  permits  phase  shifts  between  forcing  and  feedback 
processes  at  the  upper  boundary,  may  be  viewed  as  a  refined  version  of  the  classi- 

3 

cal  soil  model  described  by  Blackadar.  A  similar  model  has  been  solved  numeri- 

3.  Blackadar,  A.K.  (1979)  High  resolution  models  of  the  planetary  boundary 

layer.  Advances  in  Environmental  Science  and  Engineering,  1,  No.  1, 

J.  Pfafflin  and  E.  Ziegler,  Eds.,  Gordon  and  Breach,  50-85. 

4.  Deardorff,  J.W.  (1978)  Efficient  prediction  of  ground  surface  temperature  and 

moisture,  with  inclusion  of  a  layer  of  vegetation.  J.  Geophys.  Res., 

83  (No.  C-4);1889-1903. 

5.  Zhang,  1).,  and  Anthes,  R.A.  (1982)  A  high-resolution  model  of  the  planetary 

layer--sensitivity  tests  and  comparison  with  SESAME--79  data,  J.  Appl. 
Meteorol.,  21:1594-1609.  - 

6.  Tuccillo,  J.J.,  and  Phillips,  N.A.  (1986)  Modeling  of  Physical  Processes  in 

Nested  Grid  Model.  Technical  Procedures  Bulletin  No,  36.  3,  National 
Weather  Service,  Silver  Spring,  Md. 
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cally  in  a  multi-layer  setting  by  Cautenet  et  al.  for  a  comparison  with  experi¬ 
mental  data  in  the  African  savannah.  In  addition  to  describing  temperature  varia¬ 
tions  within  a  soil  slab  anal5rtically,  our  model  also  enables  us  to  determine  ex¬ 
plicitly  the  soil  flux  term  within  the  slab,  a  term  that  plays  a  central  role  in  the 
force-restore  method.  Furthermore,  the  perturbation  method  adopted  here  also 
points  the  way  for  linearizing  models  with  more  complicated  time-dependent  forc¬ 
ing  at  the  upper  boundary. 


2.  VIODKI.  \  riO\S 

The  variation  of  temperature  T  in  a  soil  slab  is  given  as  a  function  of  time  t 
and  depth  z  by  the  heat  equation 


a  T  _  ^  3 

3  t  3z2 


d  <  z  <  0 


where  is  the  thermal  diffusivity  of  the  soil.  We  attempt  to  solve  this  heat  equa¬ 
tion  analytically  subject  to  the  mixed  boundary  conditions: 

(i)  Energy  balance  is  maintained  at  the  upper  surface  (z  =  0), 

A  —  =  I  +  1,  -  I  -  H  -  LE;  (2) 

3  z  s  L  g  s 

(ii)  Temperature  at  the  bottom  of  the  slab  (z  =  d)  is  constant, 

T  =  T,  ,  (3) 

d 

Here,  surface  energy  balance  among  the  various  flux  terms  mentioned  in  Section  1 
is  represented  symbolically  by  Eq,  (2).  These  and  other  symbols  are  defined  in 
the  Glossary.  In  this  model,  we  have  made  the  following  assumptions: 

(1)  Tile  thermal  diffusivity  of  the  soil  ‘  is  constant  in  time  and  throughout  the 
slab. 

(2)  There  is  a  balance  of  all  lieat  fluxes  at  the  ground  surface;  temperature 
at  a  depth  d  is  maintained  at  a  constant  value  T^. 

(3)  The  initial  condition  satisfied  the  equilibrium  solution  for  tlie  given  forc¬ 
ing  at  the  boundaries. 

(4)  The  atmospheric  transmissivity  for  short  wave  radiation  may  be  approxi¬ 
mated  by  a  mean  value  i  . 


7.  t'autenct,  7  .,  Goulibaly,  V  .  ,  and  Houtin,  G.  (iriHfi)  (  alculation  of  ground 

temjierature  and  fluxes  liy  surface  models:  A  comparison  with  experimental 
fiata  in  the  African  savannah.  T ellus.  37H,  64-77. 


We  shall  assiune  further  that: 

(5)  The  solar  flux  I  in  Eq.  (2)  may  be  decomposed  into  a  daily  mean  T  and  a 
s 

diurnal  variation  1'  cosui  t.  Thus 


I  =  I  +  I '  cos  U)  t 
s 


where 


I  =  [  0)  a^i  +  a,  sin  oj  1/  tt  (Ref.  8) 

S  U  1  s 


l'  -  (sq  +  aj)  -  I 


(jj  =  cos  [  -  tan  (()  tan  6  ]  is  the  local  sunrise  hour-angle  in  radians, 
s 

(6)  Fluxes  other  than  solar  flux  may  be  parameterized  in  the  following  man- 


1,  =  e  oT  ’ 

L  a  a 

I  =  eo  T 
g  g 

H  =  C  p  C-  „u  (T  -  T  ) 
s  pa  HO  a  g  a 

LE  -  Lp  C’  u  A  fq  (T  )  -  q  I 
a  HO  a  m  ‘^s  g  ^a 

(7)  The  parameters  u  ,  T  ,  q  ,  and  A  are  externally  specified. 

We  emphasize  here  that  in  this  model: 

(1)  Phase  shifts  in  the  feedback  mechanism  are  permitted  through  terms  con¬ 
taining  Tg  in  Eq,  (2): 

(2)  The  flux  terms  I  and  I.E  are  nonlinear  functions  of  T  . 

g  g 


:i.  I.I\K\KIZ\I  ION  OK  I  NK  I  IM'KR  KOI  MIMO  rONDI  I  ION 


As  explained  above,  the  nonlinearity  in  our  inodel  is  due  to  the  iipiier  boundarv 
condition,  Eq.  (2),  Attempts  have  been  made  in  tine  past  to  linearize  the  surface 

4 

energy  equation,  I'or  example,  Deardorff  made  the  equation  linear  bv  expanding 

() 

T„  and  q  (T  )  about  a  given  instant  in  time  (time-step  n).  Pan  and  iMalirt'  expanded 
T  in  7.  and  the  potential  evaporation,  li^,  in  time  to  get  a  linear  equation.  We  shall 


8.  Iqbal,  M.  (1983)  An  Introduction  to  Solar  Kadiation.  Ac.adeniie  I'ress, 

Toronto,  Canada,  .390  pp. 

9.  Pan,  H.I..,  and  Mahrt,  E,  (1987)  Interaction  between  soil  livdrologv  and 

boundary-layer  development.  Houndary  1  .aver  Meteorol. ,  38: 1  8!'i-202. 


linearize  Eq.  (2)  by  expanding  T  and  q  about  their  respective  daily  mean  values 

6  S 

at  the  ground  surface.  Thus  we  write 


T  =  T„  +  D(z,  t). 


where  is  the  as -yet  undetermined  daily  mean  temperature  at  the  ground,  and  D 

is  the  deviation  from  T  .  Since  T  >>  D,  we  may  make  the  approximation 

S  S 


4  —  4  —4  — 

T  =  (T  +  D)^  =  T  +  4T  D. 


g  g 


g  g 


Assuming  that  the  diurnal  amplitude  of  is  1/3  that  of  T^,  we  write 

4  —4  —  ? 

+  4T  D/3, 
a  a  a 


Next,  we  expand  <4g(T)  about  <4g(T); 

+  (6) 

where  T  represents  either  T  or  T  ;  and  C  =  1  when  T  =  T  ,  and  C  =  1/3  when  T  = 

—  6  ^  S 

T^.  Here  (dq^/r*T)^  is  given  by  the  Clausius-Clapeyron  equation 

- )  =  -  .  (7) 

\st/;j;  1D:=  T“ 

and  is  given  by 

qg(T)  =  (380/p^)  exp  [17.27  (T  -  273.  16)/(T  -  35,  86)), 

where  m^  is  the  molecular  weight  of  water  vapor,  is  the  surface  pressure  and 
R*  is  the  universal  gas  constant.  Using  these  approximations,  we  may  now  write 
Eq.  (2)  as  a  linear  function  in  D: 

.  iii  =  f  +  rcosoit  +  E  o(Tf  +  4T^  D/3) 

3  cl  'J 

’  Z 


-  EO  (T +  4T^D)  -  yC  (T  -T  +  2D/3) 
g  p  g  a 


where  y  =  p  u  We  have  assumed  here  that  q(T  )  =  0.  5  q„(T  ),  that  is,  the 

a  a  ii^  o  s  a 

relative  humidity  of  the  air  is  50  percent. 

1.  I.INK\HI/.i;i)  MODKI.SOl.l  HON 


With  the  help  of  Eq.  (4),  we  write  model  Eqs.  (1)  through  (3)  in  terms  of  D: 


?  D  ,  3 

d  <  z  <  0, 

(9) 

3  t  3  * 

A  ^  ^  -  F  +  Fcoswt  -  bD, 

z  =  0 

(10) 

3  z 

D  =  -  T  . 

d  g 

z  =  d. 

(11) 

where 

F  -  r -t  L  -  -  ,C-  (T^  -T)-LE.  (12a) 

1-  g  P  g  a 


b  =  4o(.  4  2vC  p/3  +  Y  LA^ 


(12b) 

(12c) 


I.E  l.A  [q  (T  )  -  q(T  )].  (12d) 

ni  s  g  ^  a 

Thus,  in  terms  of  the  new  variable  D,  the  balance  of  heat  fluxes  at  the  upper 
boundary  is  among  a  soil  flux  term,  a  time-mean  flux  term  F,  a  time-dependent 
term  I'  coswt,  and  a  feedback  term  bl).  This  upper  boundary  condition  is  linear 
in  1)  jjrovided  tliat  T^^  is  known.  It  will  be  seen  later  in  Eq.  (20b)  that  the  value  of 
T^  is  determined  by  the  balance  of  mean  fluxes.  V\  ritten  in  the  form  of  Eqs.  (12a) 
and  (12b),  the  various  pliysical  processes  contributing  to  the  forcing  and  tlie  feed¬ 
back  terms  can  easily  be  isolated.  We  shall  illustrate  the  relative  importance  of 
some  of  these  terms  with  numerical  examples. 

The  linear  model  Eq.  (9)  has  a  solution  of  the  form 

n  -  l)e^  ^  cos  (ujt  -  u  -f  (i  z)  -t  [C  q  -♦  fit  2  cos  (wt  -  a  +  B  d)]z,  ( 13) 


where  the  constants  a  ,  B  ,  E,  t'l,  and  C  n  must  satisfy  both  the  differential  equa- 


tion  and  the  boundary  conditions.  Note  that  the  time-independent  component  C^z 
is  attributable  to  the  mean  flux  term  F  in  the  boundary  condition.  To  determine 
these  constants,  we  first  substitute  the  differentiated  form  of  Eq.  (13)  into  Eq.  (9) 
to  yield 

6  =  {a)/2K  (14) 

Next,  substituting  for  D  in  the  lower  boundary  condition  Eq.  (11)  and  evaluating 
at  z  =  d,  we  have 

C'pd  +  I)  (Cjd  4  e^'^)  cos  (wt  -  n  ^  Bd)  -  (la) 

Equating  coefficients  for  tlic  time-flc[)endent  and  time-independent  parts,  respec¬ 
tively,  we  get 


(T|-Tp/d, 


\|1  t  '  il  i. 


I'inally,  we  sul>stitute  for  H  an-l  its  derivatives  in  lap  (10)  at  the  upper  boundarv 
7.  -  0  to  obtain 

'«  tan  ^  |S  /(b  S,, )],  (U 


I)  •  I  "/[  (b  ■  S.j)  cos  I  ■  S  siti  -  ], 


(2  0a) 


where  ^  {,■  •  C'^  sin  ;  fl),  ( p  •“  cos  lid).  Thus,  liq.  (13)  is  a  solution 

to  Eq.  (9)  provided  that  Eq.  (20a)  is  satisfied  and  tlie  constants  are  given  bv 
Eqs.  (14)  through  (19).  \\  ith  the  help  of  Eqs.  (12a)  and  (16),  Eq.  (20a)  niav  be 

easily  identified  as  tlie  daily  meaTi  surface  energy-balance  equation 


>  (T  ,  -  T  )/fi  1  -  I, 
d  g  1 . 


1  T  ^  -  y  (  (T  -  T  )  -  l.i;. 
g  p  g  a 


(20b) 


Thus.  T  is  deterniined  bv  the  balance  of  flaily  mean  fluxes  at  tlie  ground  surface. 

_  «  - 

Once  T  ^  is  known,  T  (z,  t)  is  given  by  liq.  (4),  In  the  numerical  examples  given 

in  the  next  st'ction,  we  obtained  approximate  solutions  to  Eq.  (20b)  bv  first  com¬ 
puting  q^  (T  ^)  in  EE  with  a  guess  value  for  T  ^  and  rearranging  Eq.  (20b)  into  a 


quartic  equation  of  the  form 


aT^+bT  +c  =  0, 
g  g 

where 

a  =  eo  , 

b  =  yCp  -  A/d. 

c  =  XT  ./d  -  f-e  oT'‘^-yCT+  LE". 
d  a  a  p  a 

The  solution  was  obtained  by  using  IMSL  routines  and  then  choosing  the  appropri¬ 
ate  root. 

At  the  ground  surface  where  z  =  0,  Eq.  (13)  reduces  to 

D  =  D  cos  (wt  -  Cl).  (21) 

g 

Therefore,  the  temperature  at  the  ground  surface  is  given  by 

T  =  T  +  6  cos  (ojt  -  a).  (22) 

g  g 

We  now  see  that  D  is  the  diurnal  amplitude  and  a  is  the  phase  lag  at  the  ground 

surface. 

To  delineate  the  effect  of  phase  shifts  in  the  feedback  terms  in  the  upper 
boundary  condition,  we  give  here  a  solution  for  a  model  in  which  such  phase  shifts 
are  not  permitted.  In  this  case,  the  upper  boundary  condition  |Eq.  (10)]  becomes 

A  =  I  +  I'coscot  -  bD  cosoot.  (10a) 

3  z 

With  this  change  in  the  upper  boundary  condition,  the  only  changes  in  our  earlier 
model  solution  are  in  Eqs.  (18)  and  (19),  which  become 

a  =  tan  ^  (Sj/S^],  (18a) 

I)  =  l"/lb  -t  Sj  cos  u  -i  Sj  sin  a].  (IPa) 

On  comparing  these  equations  with  Eqs.  (18)  and  (19),  it  is  apparent  that  phase 
shifts  in  the  feedback  terms  tend  to  reduce  the  phase  lag  but  increase  the  ampli¬ 
tude  of  the  diurnal  changes.  Since,  as  we  shall  see  later,  3  '  ^  the  phase  lag 

given  in  liq.  (18a)  is  about  45°  . 


.“  •*•*!<*,■  «l»  ■I'.'tf.  '«f ^ 


5.  M  MERICAL  EXAMPLES 


To  get  a  feel  for  the  model  solution  of  Eq.  (13),  we  present  here  numerical 
examples  for  two  typical  cases,  a  summer  case  and  a  winter  case,  both  at  $  =3  5° 
Input  values  for  the  physical  constants  are  given  in  Table  1;  values  of  the  modeling 
parameters  for  the  two  cases  are  given  in  Tables  2a  and  2b.  Solutions  are  plotted 
as  functions  of  time  for  several  values  of  z  in  Figures  1  and  2.  For  the  summer 
case  in  Figure  1,  we  see  that  at  the  ground  surface  the  amplitude  for  the  diurnal 
temperature  wave  is  about  12.  5K,  with  the  maximum  lagging  local  noon  by  about 
an  hour.  The  amplitude  decreases  and  the  phase  lag  increases  with  depth.  The 
wave  is  completely  damped  at  the  bottom  of  the  slab  to  a  constant  value  of  about 
295K,  which  is  dictated  by  the  lower  boundary  condition.  On  solving  Eq.  (20b), 
we  found  a  value  of  300 K  for  T^,  It  is  interesting  to  see  that  at  a  depth  of  about 
60  cm,  the  daily  temperature  minimum  occurs  at  the  time  of  the  surface  temper¬ 
ature  maximum.  For  the  winter  case,  variations  with  local  time  and  depth  are 
given  in  Figure  2,  Compared  with  the  summer  case,  the  amplitudes  are  smaller 
and  the  phase  lags  are  slightly  larger.  At  the  surface,  we  have  T^  =  2771s,  T)  = 
9.5K  and  a  phase  lag  of  about  1,2  hr. 

To  check  the  realism  of  our  model  solution,  we  reproduced  in  Figure  3  some 
observed  data  cited  by  Kuo.  Comparison  of  f  igures  1  and  3  indicates  that  this 
model  is  capable  of  simulating  the  salient  features  of  the  observed  variations  dur¬ 
ing  the  day  both  in  amplitude  and  in  phase.  At  night,  the  amplitudes  are  some- 


Table  1.  Input  Numerical  t  onstants 


S  =  1 3  7  0  V\  I  n 

e  1. 0  X  10’® 

2  -  I 

111  see 

T  0.9 

-  2.2  .1  sec’’ 

K-’ 

A  -0,2 

C  =  1006  .1  kg’’ 
1) 

K-’ 

L  -  0. 95 

1,  =  4187(r)!'7.  8 

‘  0. 566  T 

.1  kg'’ 

r.  =0.6  (clear  sky) 
a  ' 

p,,  -  101. 3  kPa 

P  =  86400  sec 

IP-  287  .1 

'guess  values,  used  in 

Eq.  (12d)  only 

10.  Kuo,  II.  E.  (1968)  The  thermal  interaction  between  the  atmosphere  and  flu 
earth  and  propagation  of  diurnal  teiiiperature  ua\e.s.  .1.  Atmos,  Sci.. 
25:682-706. 


Table  2a.  Input  Values  for  Model  Parameters 


Summer 

Winter 

4)  (degrees) 

35 

35 

6  (degrees) 

21.  5 

-21.  5 

d  (meters) 

0.  8 

0.  8 

Td  (K) 

295 

•  274 

T  (K) 

299 

278 

a 

T  (K) 

301=:= 

277* 

g 

u  (ms  ) 

3 

4 

a 

(kg  kg‘^) 

0.  5  q  (f  „) 
s  g 

0.  5  q  (T  ) 
s  g 

HO 

5.  0  X  10’^ 

5.  0  X  10‘2 

A 

0.2 

0.  2 

n\ 

qg<Xg)  (kg  kg"^) 

0.  023 

0.  005 

guess  values,  used  in  Eq.  (12d)  only 

Table  2b.  Calculated  Constants  for  Summer  and 
Winter  Cases 


Summer 

Uq  (\\  111  “) 

207 

-207 

Uj  (W  ni  ■') 

752 

752 

u'  (rad) 
s 

1 .  B5 

1 .  29 

T  (W  m'“) 

352 

145 

['(W  m"- 

607 

400 

Cy  n<  m"^) 

13.  6 

5.  53 

c  J  (ill  ' ) 

0.  01 

0.  01 

ti  (ill  ^ ) 

6.  03 

6.  03 

,jq^/.iT  (T"’) 

1. 71  X  lO"’ 

4,35  X  10"“ 

SUMMER 


Figure  1,  Model  Diurnal  Temperature  Deviation  (From  the  Surface 
Daily  Mean)  at  Various  Depths  (Summer  Case).  Solar  midnight  is 
represented  as  0.  00  local  time 


what  in  error,  due  undoubtedly  to  the  fact  that  in  our  model,  forcing  at  the  upper 
boundary  is  constrained  to  a  single -component  cosine  wave, 

A  chief  merit  of  linear  models  is  the  extreme  ease  with  which  one  may  gain 
insight  into  various  physical  processes.  For  example,  the  relative  magnitude  of 
each  term  in  the  mean  surface-energy-balance  equation,  Eq,  (20b),  may  be  used 
to  assess  the  relative  importance  of  various  heat  transfer  processes  at  the  ground. 
Such  a  comparison  is  given  in  Table  3.  Here  the  dominant  fluxes  are  seen  to  be, 
as  is  well  known,  radiative  fluxes  due  to  the  sun,  the  atmosphere,  and  the  ground. 
For  typical  summer  conditions,  daily  mean  heat  transfer  due  to  nonradiational 
processes  is  small,  although  for  the  winter  case  where  T^>  T^,  sensible  heat  is 
transferred  from  the  atmosphere  to  the  ground. 

As  another  example,  we  examine  the  roles  of  the  time-dependent  terms  in  the 
upper  boundary  condition,  Eq.  (10).  This  can  be  done  by  examining  Eqs.  (18)  and 
(19).  Since  Cj  as  shown  in  Table  2b  is  typically  three  orders  of  magnitude  smaller 
than  b  ,  we  may  write  Eqs.  (18)  and  (19)  approximately  as 


WINTER 


Figure  2.  Model  Diurnal  Temperature  Deviation  (I'rom  the  Surface 
Daily  Mean)  at  X'arious  Depths  (Winter  Case).  Solar  midnight  is 
represented  as  0.00  local  time 

D  =  l'/[(b  -*  A  i:;  )  cos  i  Ap  sin  'i]-  (19b) 

Recalling  that  the  parameter  b  represents  effects  associated  witli  tlic  transient 
transfer  of  t)ieri7ial  radiation,  sensible  lieat  and  latent  lieat,  we  see  tliat  the  pliase 
lag  is  determined  principally  by  tlie  ratio  betwi-en  tlie  soil  tliermal  rapacity  and 
the  lieat  transfer  assoc  iated  with  the  fcc'dback  processes--the  smaller  the  heat  ca¬ 
pacity  and  the  Tiiore  efficient  the  heat  transfer  at  the  ground,  the  smaller  the  phase 
lag.  Similarly,  the  diurnal  amplitude  is  determined  bv  the  balance  between  the 
solar  flux  1'  and  the  non-solar  heat  fluxes  at  the  ground--for  a  given  \',  the  more 
efficient  the  transfer,  the  smaller  the  amplitude. 

A  comparison  of  the  magnitude  of  the  feedback  tc'rms  in  F.qs.  dllb)  and  (19b) 
is  given  in  Table  4.  |See  Fq.  (I2b)  for  dcTinition  of  b.  ]  1  or  summer  conditions, 
the  amplitude  and  phase  lag  are  affected  mainlv  bv  latent  and  sensible  transient 
heat  fluxes;  for  typical  winter  conditiotis  when  '  •  is  small,  thev  are  doini- 

nated  by  setisible  heat  transfer.  Note  from  Ftp  n21A  that  since  the  t ransient  sen¬ 
sible  heat  flux  depends  on  i  .  a”a^  il(>  '*  b'*'ction  of  Tj^  and  is  thus 

affected  in  this  model  more  bv  u  than  by  the  thermal  structure  i!i  the  lioundarv 

a 

layer.  The  net  transfer  associated  with  thermal  radiation,  although  a  cubic  func- 
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4 

Sunrise 


SieitsI 


Tune  (loco)) 


Figure  3.  Measured  Average  Diurnal  Soil 
Temperature  Variations  at  Various  Depths, 
O'Neil,  Nebraska  (After  Kuo**^).  Depths 
are  in  cm 


Table  3.  Relative  Magnitude  of  Various 
C  omponents  in  the  C  lear-Sky  Mean  Sur¬ 
face  Energy  Balance  Equation  (W  tti'") 


Summer 

\\  inter 

I 

353 

145 

273 

303 

-7 

-I3<i 

-310 

ti 

-IT 

s 

-3  5.  I 

1C.  7 

-T7I 

-  M  5 

-.'^6.  7 

-n 

S(Ul 

-M.fl 

-f).  3 

Table  4.  Coefficients  for  the  Various 
Feedback  Terms  in  the  Surface  Energy 
Balance  Equation  (W  m“^  K"^) 


Feedback  Term 

Summer 

V\’ inter 

longwave  rad 

4.65 

3,  62 

sensible  heat 

11,9 

17.  0 

latent  heat 

17,2 

5.  78 

Ah 

1:5.  :i 

13.  3 

tion  of  T,  is  relatively  small.  We  emphasize  here  that  whether  there  is  a  transient 
heat  gain  or  heat  loss  at  the  ground  depends  on  the  sign  of  D.  I'or  1)  >  0,  heat  is 
transferred  from  the  ground  surface  to  the  atmosphere;  for  I)  •  0,  the  opposite  is 
true.  Thus,  the  net  effect  of  the  feedback  is  negative.  That  is,  the  transfer  proc¬ 
esses  tend  to  smooth  out  the  diurnal  variations  ather  than  to  amplify  them. 

Thus,  from  an  analysis  of  the  feedback  mechanism,  we  may  identify  the  con¬ 
tribution  of  various  physical  processes  to  diurnal  temperature  variations  within  a 
soil  slab.  The  larger  tlie  value  of  T^,  the  more  important  is  the  effect  of  latent 
heat  transfer,  provided  that  u  is  not  zero  and  that  a  sufficient  inoisture  supplv 

3 

exists.  At  low  values  of  latent  lieat  transfer  becomes  less  important. 

To  isolate  the  effect  of  latent  heat  transfer,  w'e  recomputed  numerical  values 
for  Eqs.  (l.'l)  and  (20a)  for  tlic  two  cases  given  in  Table  2,  only  this  time  omitting 
terms  representing  latent  heat.  Sample  solutions  are  giyen  graphically  in  f  igures 
4  and  5.  As  expected,  botli  phase  lags  and  amplitudes  are  larger  tlian  tlie  corre¬ 
sponding  cases  shown  in  figures  1  and  2.  As  a  comjiarison,  values  of  i  and  1*)  arc 
listed  in  Table  a.  We  see  that  the  inclusion  of  latent  heat  transfer  is  important 
for  the  summer  case,  reducing  tlie  surface  mean  temjierature  bv  6K,  the  phase 
lag  by  almost  an  hour  and  tlie  amplitude  bv  6K.  In  tlie  winter  case,  its  effect  is 
less  important. 

IHSri  SSION 

Applying  the  perturbation  method,  we  have  derived  an  algebraic  expression 
for  the  time  variation  of  tenifierature  within  a  model  soil  slab  that  is  resting  on  a 
substrate  of  constant  temperature  and  is  heated  from  above  at  the  upper  surface. 
Evaluated  at  the  upper  surface,  the  expression  gives  diurnal  temperature  varia¬ 
tions  at  the  Earth's  surface.  This  model  tliffers  from  the  classical  model  given  bv 

3 

Blackadar'  in  that  phase  shift  betw'cen  the  heating  and  the  feedback  terms  in  the 
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SUMMER  LE=0. 


I'igure  4,  Model  Diurnal  Temperature  Deviation  (From  the  Surface 
Daily  Mean)  at  Various  Depths  (Summer  Case)  but  With  No  Evapora- 
tional  Cooling  at  Ground  Surface 


upper  boundary  condition  is  permitted.  As  a  result,  the  phase  angle  of  the  soil 
heat  flux  is  allowed  to  differ  from  that  of  the  solar  heating.  In  fact,  Eq.  (10) 
states  that  the  soil  heat  flux  at  the  ground  is  determined  by  three  components: 

(Da  time-invariant  component  represented  by  F;  (2)  the  diurnal  component  of  so¬ 
lar  heating,  I  'cosutt;  and  (3)  the  linearized  feedback  bD  due  to  the  diurnal  variation 
of  ground  and  air  temperatures.  Note  that  the  last  two  terms  are  not  constrained 
to  be  in  pliase  with  each  other. 

Although  phase  shift  in  the  feedback  terms  is  permitted,  the  only  external 
periodic  forcing  in  our  model  is  still  limited  to  the  I'  cos  ..  t  term.  As  a  result. 


our  solution  for  is  also  a  function  only  of  a  single  component  cosine  wave  in  t. 


llowe'ver,  if  an  external  forcing  witli  periods  wliich  are  sub-7Tiult iples  of  jieriod  I’ 
is  imposed  on  the  upper  boundary,  the  model  solution  will  then  possess  composite 


waves.  For  example,  if  terms  involving  I 
resented  by  a  Fourier  cosine  series: 

K 
I 


F' 


and  q  in  Fq. 


(2  )  could 


I'll 


I. 


,  C  T 
p  a 


Z!  f.  eos  ku)t, 
m  a  p  k 


(2  2) 
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IS 
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1 


'■S' 


£ 


It  ' 


I 

•V 


‘V 

•'■S' 


••V 

V.v 


’.V 


■sJS 


I 

-S’ 


•:-s 


<M  J 


Summer 


Winter 


Evap  /  0 

IJvap  =  0 

Kvap  /  0 

Kvap  =  0 

T  (K) 

•JOO.  0 

306.  0 

277.  0 

278.  0 

b  (W  ni'“K'^) 

:V3.  7 

16.  8 

26.4 

20.  7 

!i  ((k'f,') 

15.  8 

23.  8 

18.  5 

21.4 

I)  (K) 

12.4 

18.4 

9.  54 

11.0 

then  the  upper  boundary  condition  for  our  model  would  become 
K 

I  i  I )  ^ 

^  =  I'  n  COS  kuit  +  I  cosu)t  -  br,D, 

,z  ^  k  =  0  ^ 


gwiNBwwwpiyl 


(24) 


where 


Fq  =  I  -  C»T^^  -  vCpT^  -  vLA^q^fTgl.  (25) 

bp  =  4eoT^  +  YCp  +  YLA^(3qg/3T);jrg.  (26) 

f,  =-  /^  (I.  +  yC  T  +Y  LA  q  )  cos  ktot  dt,  k  =  0. . .  .,  K.  (27) 


In  Eq,  (27),  C  =  1  if  k  =  0  or  K;  C  =  2  otherwise.  In  practice,  if  this  soil 
model  is  coupled  with  an  NWP  model,  the  period  P  in  Eq.  (27)  could  be  defined  as 
the  24 -hr  period  preceding  the  current  time-step.  The  important  point  here  is 
that  Eq.  (24)  is  linear  in  D  and,  therefore,  Eq.  (9)  may  be  solved  in  terms  of  a 
cosine  series.  In  fact,  the  solution  is  simply  the  linear  sum  of  the  K  +  1  independ¬ 
ent  solutions,  each  satisfying  one  of  the  K  +  1  components  in  Eq.  (24); 


K 

D  =  Lj^exp  (B^z)  cos  (kwt  -  +  Sj^z) 

k  =  l 


K 

(^0  ^  ^  "  “k  ^  6j^d)]z. 


(28) 


Here  the  constants  are  determined  by  imposing  Eqs.  (24)  and  (11),  and  then  equat¬ 
ing  coefficients  for  each  of  the  components.  In  our  model,  we  have  ignored  short 
term  fluctuations  of  and  associated  with  changes  in  cloudiness.  Such  a  sim¬ 
plification  enables  us  to  retain  only  the  k  =  1  component  in  Eq,  (28),  which  is 
given  as  Eq.  (I'l), 

I  inally,  it  should  be  noted  that  this  model  may  also  be  used  to  provide  an  ana¬ 
lytical  expression  for  the  soil  heat  flux  term  used  in  the  force-restore  method.  In 
this  nietliod,  tlie  prediction  equation  for  the  temperature  of  a  slab  with  volumetric 
lieat  capacity  (  ^  is  given  Viv 


II  -  i,i; 

S 


may  thus  replace  (  'T/ .  in  (29)  with  (  •\)/ ■  7,)^j  obtained  from  Eq.  (14). 
Essentially  all  we  havt'  done  is  to  approximate  the  dce]i  soil  heat  flux  in  the  soil- 
slab  model  by  that  uiven  in  a  linearized  soil-slab  model. 

In  summary,  we  have  fornnulated  an  analytical  soil-slab  model  that  mav  be 
used  to  tmovide  a))proximatr  values  for  the  diurnal  variation  of  ground  surface 
temperature  and  to  determine  numerical  values  for  the  soil  heat  flux  term  in  the 
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force-restore  method.  Compared  with  the  classical  model  given  by  Blackadar, 
ours  has  the  following  additional  features: 

(1)  Phase  shift  is  allowed  between  forcing  and  feedback  terms, 

(2)  An  explicit  expression  is  given  for  the  diurnal  amplitude  of  the  ground  sur¬ 
face  temperature.  Although  given  in  this  model  is  still  represented  only  by  a 
single -component  cosine  wave  in  t,  it  may  be  represented  by  a  composite  of  waves 
in  cases  where  the  atmospheric  forcing  1,  ,  T  ,  and  q  are  known  functions  of  time 
for  the  entire  period  P. 

(3)  Since  soil  temperature  is  given  explicitly,  an  analytical  expression  for  the 
soil  heat  flux  term  may  be  obtained  simply  by  differentiating  Eq.  (13)  and  then 
evaluating  at  z  =  d. 

As  it  stands,  this  model  also  permits  us  to  isolate  effects  due  to  various  physical 
processes,  each  represented  by  a  different  term  in  Eqs,  (10),  (12a),  and  (12b). 

This  feature  enables  us  to  conduct  numerical  sensitivity  studies  with  ease.  For 
example,  the  response  of  diurnal  amplitude  and  phase  shift  to  moisture  availability, 
surface  roughness,  and  thermal  capacity  can  be  studied  easily  by  changing  the  val¬ 
ues  of  the  coefficients  of  individual  terms  in  Eqs.  (10),  (12a),  and  (12b). 
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temperature  at 
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